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The  1973  Yom  Kipper  War  emphasized  that  adequate  field  training 
was  essential  for  an  effective  modern  mechanized  army.  There  evolved, 
therefore,  the  requirement  (1)  to  fire  the  .50  caliber  M2HB  and  M85  machine 
guns  during  training  exercises.  As  a  result,  blank  ammunition  for  .50 
cal  weapons  is  currently  being  developed  as  part  of  the  Multiple  Integrated 
Laser  Engagement  System  (MILES)  program.  This  paper  describes  a  computer 
model  which  would  aid  in  the  design  of  the  ammunition  and  Blank  Firing 
Attachment  (BFA)  for  these  systems.  *.% 
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EXPERIMENTAL  PROCEDURE 


Weapon  Cycling 

The  .50  cal  M85  machine  gun  was  chosen  as  an  example.  Figure  1 
depicts  the  functioning  of  blank  ammunition  in  this  weapon,  which  is 
recoil  operated.  BFA 
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Figure  1.  Bolt-Barrel  Recoil. 
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Burning  propellant  in  the  chamber  produces  a  pressure  Pc.  At 
some  point  in  time,  the  mouth  of  the  cartridge  case  opens  and  gas  flows 
through  the  barrel  producing  pressure  in  the  BFA.  Chamber  pressure  in 
the  breech  and  BFA  pressure  on  the  muzzle  of  the  barrel  accelerate  the  bolt 
and  barrel. 

The  bolt  and  barrel  initially  act  as  a  single  body  during  cycling 
of  the  weapon.  But  within  a  few  milliseconds  an  accelerator  separates  the 
bolt  from  the  barrel,  and  the  weapon  cycles  independent  of  the  gas  pressure 
in  the  chamber  and  BFA.  This  model  is  concerned  only  with  interior  ballis¬ 
tics  during  the  time  the  bolt  and  barrel  remain  in  contact,  where  is  the 
effective  spring  constant  and  Ci  is  the  spring  damping  constant. 

Computer  time-displacement  for  the  bolt  and  barrel  of  the  M85 
machine  gun  when  ball  ammunition  is  fired  (2)  has  been  obtained.  The 
bolt-barrel  contact  time  lasts  only  a  few  milliseconds. 


Ball  Ammunition  Performance 


Pressure-time  traces  .50  cal  H33  ball  cartridge  were  taken  at  the 
chamber,  case  mouth  and  middle  of  the  barrel  locations.  Peak  force  due  to 
chamber  pressure  was  about  10,000  lbs;  ballistic  cycle  lasted  about  2  to 
3  milliseconds;  and  impulse  was  10-11  lb-sec.  The  blank  ammunition  and 
BFA  should  be  capable  of  duplicating  this  performance. 


Test  Set-Up 


Because  it  was  too  difficult  to  modify  the  receiver  in  the  M85 
machine  gun,  a  heavy-walled  .50  cal  test  barrel  was  modified  instead,  and 
BFA  simulator  constructed  to  obtain  pressure  traces  of  .50  cal  blank  ammuni¬ 
tion.  Pressure  stations  were  drilled  at  mid-chamber  (Pc),  case  mouth  (Pcm) * 
midbarrel  (Pmb),  and  the  BFA  (Pb)  (Figure  2). 
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Figure  2.  Cal  .50  Test  Barrel  and  BFA. 
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Shock  Waves 

Examination  of  the  pressure- time  traces  for  the  chamber,  mid¬ 
barrel,  and  BFA  (Figures  3,  4,  and  5)  shoved  shock  waves  occurring 
between  the  chamber  and  BFA, 


Figure  3,  Chamber  Pressure  Versus  Time. 


Figure  4.  Nidbarrel  Pressure  Versus  Time. 
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Figure  5,  BFA  Pressure  Versus  Time, 


Further  examination  of  the  first  return  shock  wave,  however, 
indicates  it  does  not  occur  immediately  upon  arrival  of  the  leading  shock 
at  the  BFA,  but  about  1  millisecond  afterwards  (Figure  6).  Meanwhile, 
shock  wave  theory  states  that  a  shock  wave  forms  at  the  muzzle  as 
soon  as 


Zb  ^  2  1  M*2  -  1  ±  1 

Pe  /  Y  +  1 


“  rcrit. 


(1) 


where  Me  is  the  Mach  number  of  the  gas  at  the  muzzle  and  y  is  its  ratio 
of  specific  heat. 
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Figure  6.  Shockwave  Location  Versus  Time. 

This  computer  program  monitors  the  ratio  Pfc/Pe,  and  when  it  be¬ 
comes  larger  than  rcr^t  the  gas  flow  to  the  BF A  is  cut  off.  The  pressure 
ratio  across  the  leading  shock  wave  depends  on  the  ratio  of  chamber 

pressure  to  ambient  pressure  at  the  case-mouth  opening,  on  the  ratio  of  the 
speed  of  sound  in  the  gas  to  the  speed  of  sound  in  the  air,  and  on  the  values 
of  y  in  both  gases  (propellant  gas  and  air). 


THEORETICAL  MODEL 
Three  Phases 

Propellant  combustion  and  gas  flow  occur  in  three  separate  phases 
(Figure  7).  Phase  I  occurs  in  the  chamber  and  includes: 

a.  Isochoric  or  constant  volume  combustion  until  the  mouth 
of  the  cartridge  case  opens. 

b.  Quasisteady  combustion  and  two-phase  flow. 

c.  Flame  quenching* 

d.  Quasisteady  isentropic  flow  through  a  converging  nozzle. 
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For  Phase  I  a  lumped  ballistic  model  was  used  wherein  all  gas 
properties  within  the  cartridge  are  assumed  to  be  uniform.  This  phase  also 
supplies  the  initial  and  boundary  conditions  for  Phase  II.  Phase  II  des¬ 
cribes  the  nonsteady  flow  in  the  barrel  and  includes  the  leading  shock  wave 
contact  surface  and  any  compression  or  rarefaction  waves.  The  method  of 
characteristics  is  used  to  solve  the  fluid  flow  equations  for  Phase  II. 
Finally,  Phase  III  defines  the  nonsteady  flow  entering  the  BFA  the 

pressure  buildup  in  the  BFA,  and  the  motion  of  the  bolt-barrel  assembly. 
Interpolating  along  the  characteristics  at  the  muzzle,  it  Is  possible  to 
determine  the  state  of  the  gas  entering  the  BFA.  This  flow  is  assumed  to 
cease  when  the  return  shock  forms. 


Development  of  Equations 


Phase  I.  An  Important  parameter  in  the  Phase  I  equations  is  the  solid  mass 
fraction  e  For  low  values  of  loading  density  (characteristic  of 

blank  ammunition)  the  solid  mass  fraction  c  may  be  approximated  by 


E 


£* _ 

Pp  +  Pg 


(2) 
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Accordingly,  it  is  then  possible  to  define  a  discharge  coefficient  C(j 
for  the  mixture.  _  m 

Cd  "  APc  (3) 

For  a  diabatic  nozzel  (constant  temperature) 

c  -  ■  X- -  -  .  (A) 

Y  F(Tc/Tq)  (1-e)  e 

If  one  assumes  an  initial  Kas  velocity  (v0)  due  to  gas  flow  from  the  primer 


2F(TC/T0)  (1-e) 


r  2  cT 

ll/2| 

[2  CT  -  2  (1  -  e)  F/Tfil 

Ct  Tq 

( 1  __  c\r' 

1(1  -  e)  F  IV 
^  Tq 

2  Cj  -  (1  -  e)  F/T(| 

c  CT  -  (1  -  e)  F/T0 

Vi  “ 

For  no  heat  transfer  between  the  phases 
T  o  1  1/2  p 


2  Y  T 

_  ( 

(Y  -  1)  (1 


- 1  [  LUj  :v_i-y 1 1 

-  e)  F  Tc  J  [  V  +  1/  \  Y  +  1/ 


(6) 

I  1/2 


With  these  two  parameters  defined,  it  is  now  possible  to  formulate 
a  set  of  13  differential  equations  describing  the  interior  ballistics 
(7,  8,  9)  for  Phase  I  during  the  quasisteady  combustion  and  two-phase  flow. 

A  nondeterred,  rolled-ball  propellant  is  used  in  the  propellant  charge. 

1.  Burning  rate  equation  is 

f.BP" 

dt  c  (8) 

2.  Rate  of  change  of  surface  area  for  a  propellant  grain  is 


derived  by 


-  2tt  I  tt(R  -  r)  +  4  (r  -  6) 


1  < 
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3.  Rate  of  change  of  grain  volume  equals 


dVi 
dt 

becomes 


—  -  -  2  7T  (R  -  r2)  2  -  2ir2  (r  -  6)  (R  -  r)  -  4ir  (r  -  6)2  BP®  (10) 

4.  Rate  of  change  of  number  of  propellant  grains  in  cartridge 

d  Nh  _  -e  Ch  Pc  Ar 
dt  ps  Vi 

5.  Rate  of  discharge  of  propellant  gas  from  cartridge  case  is 
d  Cc 


=  (1  -  e)  Pc  Ac  CA 


(ID 


dt  ~  -/  *C  ftc  v,d 

6*  Rate  of  increase  of  propellant  gas  in  cartridge  is  solved  by 


=  Ps  s±  Nb  BPn 


(12) 


7.  Rate  of  change  in  amount  of  gas  in  cartridge  is 


d  Cc 
dt 


d  C+ 
dt 


d  Cp 
dt 


(13) 


8.  Rate  of  change  of  gas  temperature  in  cartridge  is  derived  by 
d  Cr 


AIc.-L  [ 

dt  Cc  l 


Y  T, 


c  dt 


d  Cj 
L°  dt 


d  Cc 
lc  ~dt 


(14) 


9.  Rate  of  change  of  gas  pressure  in  the  cartridge  is 


Alt 

dt 


— — — 2  V 


c  lc 


Ale.  +  v-  C  d  Tc  -  C  T  d  Vc  I  + 
dt  +  Vc  Cc  dt  Cc  Tc  dt  1  + 


Fig_Cig 


l 


d  T, 


d  Vc 


(15) 


Tig  l  Vc  dt  Vc‘  dt 

10.  Rate  of  change  of  gas  density  in  cartridge  case  is  found  by 

d  Tc 


A !&,Ia  J_  Alfe.Ptv 

dt  F  Tr  dt  T c7 


dt 


(16) 


11.  Rate  of  change  of  solid  mass  density  is 
d£p  .  ~  Qc  Pc  Ac  ~  Cc 

dt  U0  v ’ 
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12,  Rate  of  change  of  solid  mass  fraction  is  determined  by 


dt  (pp  +  pg)  dt  (pg  +  pp)2  dt 

13.  Rate  change  of  free  volume  is 

j  vr.  -  “  iifi  dpp  _  n  d  C.c  (iq\ 

dt  ps  dt  n  dt 

These  13  simultaneous  differential  equations  were  solved  using  a  fourth 
order  Runge-Kutta  integration. 

Flame  Quenching:  Flame  quenching  occurs  when  the  pressure  starts 

to  drop  rapidly.  Subsequently,  the  flow  is  assumed  to  be  quasisteady- 
insentropic  as  through  a  converging  nozzle.  Solution  of  these  equations 
is  shown  in  Figure  8  for  the  case  where  e  *  variable  (vQ  ^0)  and  z  =  0.5 
(vQ  0).  Figure  8  shows  the  experimental  curve.  Since  a  BFA  was  not  used, 
no  shock  waves  occurred. 
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Gas  Leakage.  Gas  leakage  in  the  BFA  has  a  significant  effect  on  the  pressure. 
This  leakage  (Figure  9)  occurs  through  the  forward  end  of  the  BFA  (ABFA) , 
through  the  clearance  between  the  BFA  and  the  barrel  (^BLK)  and  after  the 
barrel  has  recoiled  sufficiently  through  the  vent  port  ^PRT. 


Figure  9.  Cal  .50  Barrel  and  BFA. 

Equation  of  Motion.  The  equation  of  motion  for  the  bolt-barrel  assembly  is 
determined  from  applied  forces.  Pressure  in  the  BFA  and  in  the  cartridge 
case  supply  the  accelerating  forces.  Five  simultaneous  differential  equa¬ 
tions  must,  however,  be  solved  to  determine  this  motion. 

1.  Equation  of  state  for  the  BFA  is  solved  by 

d  _  mbR _  d  Th  .  R  Th  dmh  RTb  mb  d£  ,2(), 

dt  Ub  +  Adlib  dt  Ub  +  AdUb  dt  "  (U^  +  Ad^)2  dt 


2.  Energy  equation  for  BFA  becomes 


me  (1  -  3) 


ve2 

p  m  ,  VC 

/  abfa\ 

Ph  Ab  dib  ,  d 

p  e  *  2gJ 

V  Ac  J  = 

J  dt  dt 

^7  +  77  <"bgCvTb)  + 


cp  Tb  ®ob 


or  solving  for  gas  temperature 


d  Tb  _  me  (1  -  3) 
dt  Cv  mb 


Te 


i1 


(21) 


ABFAN  Pb  Ar  dib 

Ac  J  gJCymb  dt 


Cp  Tb  mob  _  Tb  mb 
mb  Cv  b 


(22) 
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where  me  is  the  rate  at  which  gas  enters  the  BFA  and 

j  y  +  1 

mob  =  Pb  (ABLK  +  NH  APRT)  "v/r1^  (t^i)  Y  "  1  ^ 

3.  Equation  of  motion  is  found  by 

dv  dv 

Pc  Ac  +  Pb  Ab  -  Kx  (Ab  +  x0)  -  Cl  V  -  Fx  -  F2  =  Mt  ^  +  16.56  ^ 


4.  Conservation  of  propellant  gas  mass  is  computed  by 

•  •  • 

=  me  -  mob  (25) 

5.  Equation  of  barrel-bolt  displacement  is 


BFA  Pressure-Time  Curves.  The  computed  and  experimental  pressure-time 
curve  for  the  BFA  and  the  experimental  pressure- time  curves  are  shown  in 
Figure  10.  The  effect  of  the  reverse  shock  wave  is  evidenced  by  the  decrease 
in  pressure  following  the  peak.  Meanwhile,  the  rise  in  pressure  for  the 
experimental  curve  following  the  peak  is  due  to  the  arrival  of  the  return 
shock  from  the  breech. 
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Figure  10.  BFA  Pressure  Versus  Time. 
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CONCLUSIONS 

1.  This  model  appears  adequate  for  describing  the  quasisteady  combustion 
and  two  phase  flow  from  the  .50  cal  blank  cartridge.  It  can  also 
simulate  the  leading  shock  wave,  compression  wave  and  pressure  buildup 

in  the  BFA. 

2.  The  effective  time- averaged  solid  mass  fraction  e  of  about  0.5 

seems  reasonable  for  simulating  the  peak  pressure  in  the  cartridge  case. 

3.  Quenching  of  the  propellant  flame  immediately  following  the  start 

of  depressurization  of  the  cartridge  accounts  for  the  shape  of  the 
pressure- time  curve  in  the  cartridge  following  peak  pressure. 

4.  The  pressure  drop  in  the  BFA  chamber  occurs  when  the  first  return 

shock  wave  forms  and  travels  back  toward  the  breech. 

5.  The  method  of  characteristics  can  be  used  to  solve  equations 
describing  nonsteady  isentropic  gas  flow  from  cartridge  case  to  BFA. 
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LIST  OF  SYMBOLS  O) 

A^  =  cross-sectional  area  of  BFA 
Cp  =  specific  heat  for  gas 
Cg  =  specific  heat  for  particles 
CT  =  (1  -  e)  Cp  -  eC$ 

Ci  »  spring  damping  constant 
e  =2.718  —  base  of  natural  logs 
F  -  propellant  impetus 
T?l  =  bolt  receiver  friction  force 
F2  =  barrel  receiver  friction-force 

16.50  Kjj  =  effect  of  number  of  rounds  in  the  ammunition  belt 

Ki  =  effective  spring  constant  for  bolt-barrel  assembly 

&b  =  travel  of  bolt-barrel  assembly 

mt,  =  mass  of  the  gas  in  BFA 

Mt  =  total  moving  mass 

pb  =  pressure  in  BFA 

R  =  gas  constant 

Tc  =  gas  temperature  in  the  cartridge 

T0  =  isochoric  adiabatic  flame  temperature 

Ub  =  initial  volume  of  BFA 

XG  =  initial  compression  of  spring 

r)  =  gas  covolume 

pg  =  density  of  gas  in  the  mixture 

pp  =  density  of  solid  particles  in  the  mixture 
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